Dynamical lattice instability versus spin liquid state in a frustrated spin chain system 
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The low-dimensional s — 1/2 compound (NO)[Cu(N03)3] has recently been suggested to follow 
the Nersesyan-Tsvelik model of coupled spin chains. Such a system shows unbound spinon excitations 
and a resonating valence bond ground state due spin frustration. Our Raman scattering study of 
(NO)[Cu(N03)3] demonstrates phonon anomalies as well as the suppression of a broad magnetic 
scattering continuum for temperatures below a characteristic temperature, T < T* ~ 100 K. We 
interpret these effects as evidence for a dynamical interplay of spin and lattice degrees of freedom 
that might lead to a further transition into a dimerized or structurally distorted phase at lower 
temperatures. 

PACS numbers: 75.10.Kt, 75.10.Pq, 78.35.+C 
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I. INTRODUCTION 

Quantum spin systems with competing interactions 
and reduced dimensionality of the spin exchange show 
exotic physicsii In particular, there is the opportunity 
to observe unusual spin liquid ground states, exotic frac- 
tional excitations and to test various theoretical concepts. 
This has triggered intensive research activities. However, 
the realization of such states is challenging, as more con- 
ventional, long range ordered or dimerized (spin-Peierls) 
phases exist as competing instabilities. The needed fine 
tuning and tailoring of the crystallographic and electronic 
properties requires a considerable knowledge of materials 
aspects as well as deep theoretical understanding. So far, 
most of the studied systems can be classified as quasi 1- 
dimensional (ID) quantum spin chains, coupled chains, 
2-dimensional (2D) triangular and kagome lattices. 

The recently synthesized low dimensional compound 
nitrosonium nitratocuprate (NO)[Cu(N03)3] has been 
discussed in the framework of a frustrated, quasi-2D sys- 
tems without long range magnetic orderingi^ The crystal 
structure of (NO)[Cu(N03)3] consists of layers of CuOe 
octahedra that are coupled along the crystallographic 
6 direction by triangular (N03)~ (nitrate) groups, see 
the inset in Fig. [T] These chains are weakly coupled 
via double groups of singly connected nitrates. In ad- 
dition (NO)"^ (nitrosonium) cations are located between 
the chains. Interestingly, these three structural subunits 
not only lead to exotic magnetic properties but also to 
phonon modes with three frequency ranges and a charac- 
teristic low energy onset, similar to fingerprints, as shown 
further below. 

The magnetic susceptibility, specific heat, and electron 
spin resonance (ESR) measurements in (NO)[Cu(N03)3] 
evidence pronounced spin fiuctuations without any hint 
for long-range magnetic order down to 1.8 K;^ A broad 
maximum in the magnetic susceptibility x(T) is observed 
around T* ^ 100 K, characterizing the energy scale of the 
spin exchange. There is no evidence for an excitation gap 



of the triplet states (spin gap) that would lead to a sharp 
suppression of x(2^) ^t lowest temperatures. However, 
x(r) with T <T* is smaller than expected for an isolated 
s = 1/2 Heisenberg chain^i^ 

From a detailed analysis of x(T) it has been deduced 
that (NO)[Cu(N03)3] is a realization of the 2D confeder- 
ate flag or Nersesyan-Tsvelik modeli^i^ This model con- 
sists of spin chains with additional direct and diagonal 
coupling between neighboring chains. The diagonal cou- 
plings are related to spin frustration as not all exchange 
paths can be satisfied simultaneouslyii The particular ex- 
change geometry is known to lead either to a resonating 
valence bond (RVB)i or a valence bond crystal (VBC) 
ground state with a macroscopic ground state degen- 
eracy. These states support free, propagating spinous 
and more exotic topological excitations. The assign- 
ment of the compound to the Nersesyan-Tsvelik model 
is based on a critical ratio of three exchange coupling 
constants: the leading magnetic exchange coupling along 
the chains, J=170 K, the interchain coupling, J', with 
—8.5 K < J' < 15.3 K, and the exchange interaction J2 
along the diagonals. The latter is proposed to be exactly 
1/2 • J, which leads to strong spin frustration, suppresses 
long-range order (LRO), and allows quantum fluctuating 
and entangled ground states. However, in the proximity 
to the critical coupling also dimerized phases exist 1^ 

In a recent density-functional-theory (DFT) band 
structure calculation doubts about the Nersesyan-Tsvelik 
model for (NO)[Cu(N03)3] were raised^ Considering a 
different geometry of superexchange pathways, a similar, 
leading antiferromagnetic exchange coupling was found 
with J = 200 K. However, much weaker couplings with 
J' = 2 K link the chains into layers. Therefore the com- 
pound would be better described as a uniform spin-chain 
with weak and nonfrustrated interchain coupling. Irre- 
spective to the DFT the present thermodynamic experi- 
ments do not allow a clear consensus which model leads 
to a better description. 

For a deeper understanding of the title compound it 



is of pivotal importance to probe both lattice and mag- 
netic degrees of freedom using a spectroscopic experimen- 
tal techniques.— 1^ In addition, a previous Raman scat- 
tering study of (NO)[Cu(N03)3] reported only four of 
the expected 45 phonon modes which further motivates 
a refined studyj^ With our Raman experiment we sup- 
port a frustrated spin model with a coupling constant of 
Jb« 150 K and give further evidence for strong phonon 
anharmonicity in (NO)[Cu(N03)3] possibly leading to a 
dynamic dimerization. 



II. EXPERIMENTAL DETAILS 



of nitrosonium nitratocuprate 
were obtained by means of a wet 



Single crystals 

(NO)[Cu(N03)3 
chemistry routci^ The crystals are unstable under 
ambient conditions and were therefore kept in vacuum to 
prevent degradation of the freshly cleaved crystal surface 
during the investigation. Raman scattering experiments 
were performed in quasi-backscattering geometry using 
a A = 532 nm solid state laser. The laser power was 
set to 3 mW with a spot diameter of approximately 100 
/im to avoid heating effects. All measurements were 
carried out in an evacuated closed cycle cryostat in the 
temperature range from 10 K to 295 K. The spectra 
were collected via a triple spectrometer (Dilor-XY-500) 
by a liquid nitrogen cooled CCD (Horiba Jobin-Yvon, 
Spectrum One CCD-3000V). 



III. EXPERIMENTAL RESULTS 

At temperatures above 160 K, the crystal structure 
of (NO)[Cu(N03)3] is monoclinic (space group P2i/m) 
with two formula units per primitive cellf^ii^ For this 
crystal symmetry the factor group analysis yields 45 
Raman-active phonon modes: Teaman — 25 • ^g 4- 20 • Bg. 
The corresponding Raman tensors are given by: 
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Raman scattering experiments performed in different po- 
larizations at room temperature (T = 295 K) within the 
ab plane of the crystal do not reveal a noteworthy depen- 
dence on the polarization configurations. This is prob- 
ably due to multiple reflections of the incident light in 
the transparent crystals. Therefore we have removed the 
analyzer to obtain maximum scattering intensity with 
(xx + yx) polarization in all following measurements. 

Raman scattering spectra collected at 295 K and 10 K 
in the frequency range up to 2400 cm""'^ are presented in 
the left panel of Fig.[T] According to the three structural 
subunits the lattice vibrations in (NO)[Cu(N03)3] can 
be separated into distinctive groups of frequencies corre- 
sponding to external translational and rotational modes 
(70 - 200 cm-i), internal modes of CuOe (200 - 600 




1500 2250 2275 100 200 

Raman shift (cm') 

FIG. 1. (Color online) Left panel: Raman spectra of 
(NO)[Cu(N03)3] obtained in [xx) polarization at 10 K and 
295 K. The filled background emphasizes the temperature 
evolution of the observed scattering continuum. The crys- 
tallographic structure of (NO)[Cu(N03)3] is shown as an in- 
set within the a6-plane. The N(II)OJ nitrate groups connect 
Cu chains along the &-axis, while two N(I)Oj nitrate groups 
form more complex but weaker interchain connections along 
the a-axis. Right panel: Temperature evolution of the Raman 
spectra in the low-frequency range in (xu) polarization. 
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FIG. 2. (Color online) Temperature dependence of the low 
energy phonon frequencies. The red squares correspond to 
phonon modes that appear at low temperatures, T < T* ^ 
100 K, the latter given by a dashed vertical line. Lines are 
guides to the eye 



cm~^), nitrate groups NO3" (700 - 1500 cm~^)ii as well 
as NO"*" groups^ (w2280 cm~^). Fitting the spectrum 
with Lorentzian profiles allows to distinguish 41 phonon 
modes in the entire frequency region. This is in good 
agreement with the factor group analysis. It is notewor- 
thy that there exists a well defined low energy onset of 
the phonon modes at Y^iow ~70 cm~^ = 100 K, which co- 
incides with the above mentioned T* from x(^)- There is 
a pronounced temperature dependence of the phonons as 
well as of a broad scattering continuum. In the following 
we will study both effects in detail. 

In the right panel of Fig. [T] the temperature develop- 
ment of the spectra in the low frequency range is plot- 
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FIG. 3. (Color online) Analysis of the phonon frequency shifts 
(symbols) using a model (lines) based on phonon-phonon de- 
cay processes.^ The corresponding modes show deviations 
for temperatures T <T* , especially for higher frequencies, as 
described in the text. 



ted. Strong changes as function temperature involve 
frequency, intensity and linewidth. In order to analyze 
and separate the temperature dependence of the phonon 
modes from the broader magnetic scattering, all spectra 
were analyzed by fitting them to lines with Lorentzian 
profiles. 

In Fig. [21 we illustrate shifts in phonon frequency in the 
low energy region as function of temperature. Focussing 
here on the modes at around 190 and 220 cm~^ (right 
panel) it is obvious that there exists a gradual structural 
evolution. The onset of this effect is above 200 K. The 
strongest changes including the appearance of new modes 
are at T < T* w 100 K. This is the characteristic tem- 
perature where the magnetic susceptibility x(2^) shows 
a broad maxinmm. At low temperatures the total num- 
ber of phonons exceeds the expected 45 modes from the 
factor group analysis, i.e. we have to consider a reduced 
local symmetry in (NO)[Cu(N03)3]. 

Figure [21 left panel, shows the temperature evolution 
of some external translational and rotational phonon fre- 
quencies together with a fit based on phonon-phonon de- 

^uph{T)=u:o + c(l+ 2 



FIG. 4. (Color online) Analysis of the phonon linewidths 
(FWHM, symbols) using a phenomenological model (lines). 
Several modes show deviations for temperatures 200 K < T < 
T*. 
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This fit approximates the anharmonic phonon contribu- 
tion. The experimental data deviate in the low temper- 
ature regime, setting in for T <T* and the modes have 
a lower frequency than expected. Also internal modes of 
CuOg, expected in the range 200 - 600 cm~^, show an 
anomalous temperature dependence, see right panel of 
Figure [3l Noteworthy is that positive as well as negative 
deviations from the anharmonic fit exist. 

The linewidths (full width at half maximum, FWHM) 
of selected modes in the frequency region of external (in- 
ternal) vibrations are depicted in the left (right) panel 
of Fig. m Our analysis shows that many phonons are 
anomalous, i.e. they show variations as a function of 
temperature that are not consistent with anharmonicity. 
The solid lines in the temperature dependent linewidths 
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FIG. 5. (Color online) Integrated intensity over temperature 
of selected phonon modes. Lines are guides to the eye. 



are a fit using a phenomenological model: r(T) = 

^O cxp(rJo/fcBT)-i ' "^i^^ ^J ^'^i'^S a mode dependent fit pa- 
rameter. The phonon linewidths show a general decrease 
with reduced temperature. This is due to anharmonicity 
and phonon-phonon scattering. Noteworthy, however, is 
the superimposed change of curvature leading to an ad- 
ditional drop of linewidth for T < T* « 100 K of the 
modes at 140, 154 and weaker for the one at 65 cm~^. 

Fig. [5] plots the Bose corrected integrated intensity, in 
analogy to frequency and linewidth of the phonon modes. 
It is noteworthy that the intensity evolution is different 
for each mode without an obvious dependence on the fre- 
quency of the mode. The observed effect is therefore not 
related to the scattering volume, i.e. by a change of the 
optical penetration depth. It has rather to be individu- 
ally assigned to certain spin-phonon coupling parameters 
as discussed below. The integrated intensity of the inter- 
nal CuOg modes (right panel. Fig. [S]) shows evidence for 
a second characteristic temperature at 200 K marked by 
a dashed lines. 
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FIG. 6. (Color online) Magnetic Raman scattering with a fi- 
nite energy maximum EmaxiT) determined from fitting the 
data after subtracting all phonon lines. The curves were ob- 
tained at 10, 20, 40, 60, 80, 100, 125, 150, 175, 200, 225, 250, 
and 295 K. The upper inset shows EmaxiT) and the lower 
inset the integrated intensity lintiT). 



The integrated Raman scattering intensity of a phonon 
mode consists beside the scattering volume of two con- 
tributions J^ The dipole matrix elements arise from a 
phonon-induced polarization of the wave function and 
these elements vary exponentially with the interatomic 
distance. The other contribution is related to displace- 
ment induced changes in the band energies. As the effects 
in (NO)[Cu(N03)3] are specific for each phonon line they 
can not be explained in terms of lattice thermal contrac- 
tion. Therefore, dynamic changes of ionic distances are 
the most likely origin for the observed variation of phonon 
intensity. 

There exists a further Raman scattering signal in 
(NO)[Cu(N03)3] at T > r* with a much larger linewidth 
than the phonon scattering, see the solid background in 
the left panel of Fig. [TJ This continuum shows a finite 
energy maximum {Emax ~ 200 cm^^). We have no 
evidence for disorder that could broaden phonon lines. 
Therefore we attribute the continuum to magnetic Ra- 
man scattering.— 

The intensity of the continuum is strongly reduced 
at low temperatures, see lower inset of Fig. [HI We 
have to admit that the continuum overlaps with many 
sharp phonon lines. Therefore it is difficult to determine 
whether its intensity is reduced exactly to zero or to a fi- 
nite magnitude. The upper inset shows the energy of the 
maximum. There is a continuous shift by approximately 
10% of its position to higher energy with decreasing tem- 
peratures. Superimposed on this shift there is a jump at 

In a ID s = 1/2 spin system Raman scattering on 
spinon excitations may be observed as a broad maximum 
with a typical energy of Emax = 2 J, with J the exchange 



coupling along the chain direction. From the maximum 
in Fig. [6] we derive a coupling constant of J(,«150 K. This 
value compares very well with previous values from ther- 
modynamicSj^ J=170 K, and still reasonably well with 
Jb~200 K from DFT calculations^ Such magnetic scat- 
tering is only observed if the Raman exchange scattering 
Hamiltonian H^ = ^ s ■ s does not commute with the 
magnetic exchange Hamiltonian of the spin systems, i.e. 
additional couplings causing frustration or a dimerization 
is mandatory.— 

The strong decrease of scattering intensity for T < T* 
K is not expected. In contrast, spinon and two-magnon 
scattering is only weakly temperature dependent and 
only broadened by thermal fluctuations with T compara- 
ble to the coupling constant. This has been shown, e.g. 
for magnetic Raman scattering in the kagome latticej^ 
The temperature effect could be related to a suppression 
of the free spinous by opening a static or dynamic gap. In 
the spin Peierls system CuGe03 a similar depression has 
been observed with decreasing temperatures in proximity 
to the dimerized phase»ii Although x{T) is smaller than 
expected, there is no pronounced effect as expected due 
to the freezing out of triplet excitations. Nevertheless, 
the frequency shift of the continuum seams to indicate a 
continuous change of the exchange parameters with T in 
the overall temperature range as well as a discontinuous, 
first-order like change at T*. 



IV. DISCUSSION 

In the following we will discuss the implications of the 
phonon anomalies for the structure of (NO)[Cu(N03)3]. 
Then possible magnetic ground states consistent with the 
observed magnetic Raman scattering will be evaluated. 

In discussing phonon anomalies in (NO)[Cu(N03)3] it 
is important to realize that its structural data^ evidences 
a remarkable flexibility of the NO;^ and N0+ subunit 
positions, see inset of Fig. [TJ While the Cu positions 
centered in CuOe octahedra are firmly fixed and may play 
a stabilizing role, other structural units are located at the 
general 4/ (x, y, z) crystallographic positions in which all 
X, y, z coordinates are free. This concerns all atoms of the 
two N(I)0;^ subunits or N(II) and 0(22) atoms from the 
N(II)0;^ subunits at the 2e (x, 1/4, z) positions in which 
two X— and z-coordinates are free, and all atoms from 
N0+ subunits. Here we use the notation from Rei^ for 
the labeling of oxygen and nitrogen atoms. The NO3" 
complexes are therefore allowed to rotate as a whole and 
NO"'" complexes can rotate within the ac-plane without 
changing any Cu-0 bond length. 

A thermal population of the low frequency external 
translational and rotational modes of the NO;^ and N0+ 
subunits thereby enhances the anharmonicity of all op- 
tical phonons as decay channels are opened. As these 
phonon modes have a low energy onset at Ejou; ~ 70 
cm~^ = 100 K. most of the phonon anomalies are char- 
acterized by this temperature as well. Such a mechanism 



via lattice anharmonicity also takes into account that the 
observed frequency shifts are rather moderate, around or 
below the 1% range. Therefore we suggest a tendency for 
local or fluctuating lattice distortions in (NO)[Cu(N03)3] 
that could eventually lead to a crystallographic distortion 
with reduced site symmetry. 

We have evidence for such a kind of dynamical insta- 
bility also from the high frequency phonon modes. The 
existence of a few dynamically stable positions at high 
temperature is observed as a splitting of the characteris- 
tic NO"*" mode at around 2000 cm~^. This splitting dis- 
appears at low temperatures (see high frequency range 
in the left panel of Fig. [T]). Moreover, the phonon fre- 
quency 7(N0'^) is located at low temperatures at 2280 
cm~^. This frequency reflects an absence (or minimum) 
of interaction between the NO''' units and Cu atoms. In 
the opposite case, a significant decrease of the 7(N0+) 
frequency should be observed as was demonstrated for 
several NO complexesj^^ 

Taking this structural flexibility and dynamics into ac- 
count we differentiate two effects on the magnetic proper- 
ties of (NO)[Cu(N03)3]. CuOe octahedra determine with 
their Cu-0 bond lengths the strength and with their dis- 
tortions a possible dimerization of the exchange coupling 
along the &-axis. This coupling establishes the ID mag- 
netic correlations of this compound, which is dominat- 
ing. In contrast, the NO^ complexes mediate the direct 
and diagonal frustrating couplings of the spin chains with 
each other 1^ These couplings are critical for establishing 
different kinds of quantum disordered ground statesi^, 
and whether (NO)[Cu(N03)3] realizes the Nersesyan- 
Tsvelik model. 

The phase diagram of 2D 3\l Ji quantum spin systems 
has been investigated theoretically^^ as function of frus- 
trated couplings and thermal fluctuations in detail. De- 
spite several open questions it is clear that quantum dis- 
ordered phases, e.g. of the Nersesyan-Tsvelik model, have 
a rather small phase space. Such phases are enclosed by 
magnetically ordered or dimerized states with different 
low energy excitations. Thermal fluctuations destabilize 
long range order and thereby enlarge the phase space of 
the exotic phase^i^. We attribute the intensity variation 
of the spinon continuum to a temperature induced shift 
of the system to the boundary of an ordered state. In this 
sense the behavior of (NO)[Cu(N03)3] is unconventional 
and could be compared with a spin chain system that is 



coupled to phonons. In the latter system a dimerization 
transition (e.g. a spin-Peierls transition at Tsp) opens a 
spin gap and suppresses spinon excitations at low temper- 
atures. Such a gap, however, concerns both singlet and 
triplet excitations and it is at a temperature much smaller 
than the characteristic energy scale of the spin system. 
In (NO)[Cu(N03)3] we observe the onset of the suppres- 
sion of spinon scattering at around 250 K, a temperature 
of approximately 2J(,. A modeling of (NO)[Cu(N03)3] 
should therefore consider a 2D exchange topology includ- 
ing lattice fluctuations. A purely ID spin chain model is 
not in agreement with the observation of a spinon con- 
tinuum in Raman scattering. Additional spectroscopic 
and structural investigations of (NO)[Cu(N03)3] also at 
lower temperatures, e.g. using neutron scattering, would 
be needed to further clarify this point. 



V. SUMMARY 

In conclusion, we observe phonon anomalies in 
(NO)[Cu(N03)3] promoted by the weakly coupled struc- 
tural building blocks of the compound. The characteris- 
tic temperature of these anomalies and of the magnetic 
susceptibility coincide. In addition, a broad, temper- 
ature dependent scattering continuum of magnetic ori- 
gin is observed and attributed to free spinon excitations. 
The position of this maximum is related to the major 
exchange coupling Jh«150 K along the chain direction 
in good agreement with previous investigations. We at- 
tribute further anomalies to a pronounced anharmonicity 
and dynamics of NO;^ and NO"*" subunits that modu- 
late the coupling between the chains. Our observations 
do not directly contrast with the proposed Nersesyan- 
Tsvelik model. However, a refined description of the 
compound should take spin-phonon coupling and lattice 
dynamics into account. 
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